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fe LRODUCTION 


The perishable product is beginning to receive attention 
in the inventory problem. Often, it has a fixed lifetime 
after which it becomes useless to satisfy demands. News- 
papers, medical drugs, human blood, photographic films and 
some foodstuffs are common examples. 

An inventory problem for perishable products is often 
treated as a single-period inventory problem, sometimes 
called a newsboy problem. The newsboy problem was introduced 
by Morse and Kimball in 1950 and the development of it has 
been presented in several journals and publications. A 
brief history of it was given by Masuda [Ref. l]. 

A well-known characterization of the newsboy problem is: 

The newsboy must decide how many papers to purchase 
for resale during the day. If he buys too many papers-- 
that is, more papers than people will purchase from him-- 
he will have papers left at the end of the day and will 
incur a loss on each, since old newspapers have little 
salvage value. If, on the other hand, he purchases too 
few papers for resale, he will sell out early and his 
later customers will not be able to buy from him; he thus 
incurs a cost (lost profit) for each paper he could 
have sold, but didn't. The newsboy's problem is to de- 
cide how many papers to have on hand at the beginning of 
the day so as to minimize the losses (due to surplus or 
shortage) at the end of the day. [Ref. 2] 

Finding good solutions to the newsboy problem depends 
upon the amount of information about the magnitude of demand 
X. 


In recent research dealing with newsboy ordering deci- 


Slons where little information existed about demand, it was 





assumed that at the end of each period the demand for that 
period would be known [Refs. 1,3,4]. In order for this to 
happen, one would need unlimited supply or some way to 

measure the demand from unsatisfied customers. In practice, 
this is often not the case. For example, when the newsstand 
is empty, the later customers will not stop at that news- 
stand. Accordingly, we do not know how many more customers 
would like to buy that newspaper. Thus, we may have frequency 
data on demand only up to the value of supply. 

This thesis is intended to provide methods permitting 
solutions to the newsboy problem when we only have demand 
frequencies up to the supply value. 

In Chapter II, two procedures will be proposed to esti- 
mate the maximum demand. The first procedure uses the Gompertz 
Curve since the nature of the demand cumulative distribution 
function (c.d.f.) can be viewed as a growth curve with the 
market saturate as an asymptote. The second procedure sug- 
gested is the four-parameter Schmeiser-Deutsch Distribution 
since it can take on a wide variety of shapes so as to match 
many distributions. Also, the Schmeiser-Deutsch Distribution 
was found to be a suitable procedure for the case of limited 
supply. 

Using the estimation of maximum demand together with the 
unit cost of surplus and the unit cost of outage, the estima- 
tion of optimal order quantity is discussed in Chapter III 


for the "under uncertainty” and "under risk" cases. 


EO 





Conclusions and applications are given in Chapter IV. 
A review of the Schmeiser-Deutsch Distribution is presented 


in Appendix A. 
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II. ESTIMATION OF MAXIMUM DEMAND 


Care should be taken to estimate the maximum demand since 
the optimal ordering quantity for uncertainty depends on it. 
In this chapter, two alternatives are suggested to estimate 
the maximum demand when we have unlimited supply. Then, we 
will show how one of these methods may be used to estimate 
maximum demand when there is no data on unsatisfied demand 
(the case of limited supply). 

The procedures that appear suitable and realistic are 
(1) using the Gompertz Curve and (2) using the Schmeiser- 
Deutsch Distribution. We will discuss the Gompertz Curve 


porst. 


A. USING THE GOMPERTZ CURVE TO ESTIMATE THE MAXIMUM DEMAND 
In the terminology of economics, the demand c.d.f. can 

be considered as a "growth curve" since it has a "market 

saturate" as an upper asymptote. The Gompertz Curve is the 


most generally used growth curve. Its formula 


where X is the time series and Y is the sales amount, de- 
pends on the values of the parameters k, a and b (or their 
logarithms), it may take on any one of a variety of shapes 


[Ref. 5]. 
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Figure 1. The Gompertz Curve with 0 <a<1l, 0<b<l 


We are interested in using the Gompertz Curve with 0<a< 1 
maa 0 < b < 1 to approximate the c.d.f. of demand. This form 
of the Gompertz Curve is known to be asymptotic to k. When 
we fit the Gompertz Curve to the observed cumulative fre- 
quency of the data, k should theoretically be 1.0. In prac- 
tice, however, it may depart slightly from 1.0. To compensate 


for this, we use the following form to approximate the c.d.f.: 


The following description is taken from [Ref. 5]: 
The parameters k, a and b are obtained as follows: 
1) Be certain that the number of observations is divisible 
by 3, i.e., there are 3N periods of base data. 
2) Take the common logarithms of e@ach Y. 


meeaa the first N log Y's to obtain £, log Y; then the 


i 
second N data points to obtain Zo log Y; and the last N 
Beamnes to obtain Z 3 Leg Y < 


4) Set the smallest observation to be unity. 
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5) Substitute in the following formulas: 





AN _ Z,log Y - 2 1log Y a 
Z,1og Y- Z log Y 
HOgua = “(s-log Y = £.leog Y)- 2) (3) 
J 2 ere Aula 2 
(b -1) 
and 
i yal 
log k = = Gr log vr Jee lOG <a). (4) 
N a ol 


With values of k, a and b, one can compute Y for any 
value of X using the formula of the Gompertz Curve. Also, 
the formula can be derived from the Gompertz Curve to 


generate X from the Gompertz Curve as: 


eC LOCK 
a log ( log a ) aa 
log b 


In this chapter, we wish to estimate the maximum demand. 


If the demand distribution has an infinite upper tail, we 


h 


can not hope to estimate the 100° percentile of demand. 


We will be content, therefore, with the gg th 


percentile of 
demand. This is found by setting (Y/k) = .99 or Y = .99k 


in equation (5), i.e., 


eg (aoe 
= log a 
ee will be = leq 5 (6) 
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which may be considered to be the estimator of the maximum 
demand if the smallest observed X is unity. (If not, the 
maximum demand should be adjusted to be X plus the smallest 
X minus one.) 

As an example assume that the demands for newspapers 
meave been 10, 6, 9, 7, 5, 13, ll, 7, 8, 8. Using the Gompertz 
Curve to estimate the maximum demand, we can preceded as 
follows: ~ 

Let X be the demand for newspapers, and Y be the cumula- 


tive relative frequency of X. 


TABLE I 


DEMAND AND ITS CUMULATIVE RELATIVE FREQUENCY 


us 2 
5 Orel 
6 or 
7 0.4 
8 0.6 
9 C7 
10 0.8 
11 0.9 
1D 0.9 
13 ie C 


Since Y is a growth curve, it can be approximated by the 


X 
Gompertz Curve, Y = ka? . To determine the appropriate 


values of k, a and b using the method described above: 
1) From these data, there are 9 distinct values between 
the smallest observation (5) and the largest observation (13), 


mu? 1S divisible by 3. (Here, N = 9/3 = 3.) 


lies. 





Seen the third column of Table II, Y's are converted to 
eg Y's. 

3) Values of Z log Xe; ZLog 16 ENgKS, 
ime, fourth column of Table II. 


310g Y are shown in 

4) Since the smallest observed X is not 1.0, the maximum 
demand should be adjusted by adding the difference between the 
smallest observed X and 1.0. This is because the Gompertz 
Curve assumes the smallest observation is unity. 


5) From equation (2); 


oom geet 73 7) 7 ((—. 4737) -(-2.0969)) 
= 0.2388 
Therefore, = 
5 = mys 


From equation (3); 


ocean = c= .47 07—(—-2,0969) ) (.6175-1) 7(.2355-1) 
—ea 0 O22 
Therefore, 
a = 0.0867 


From equation (4); 


log k 


=(-2.0969-(.2355-1) (-1.0622) /(.675-1) ) 


OVS ROIS) FF 


INS 





Therefore, 


Koes. 0205 


Xx 
I.e., the Gompertz Curve is Y = pumaemennqomosieg) ‘°° °+/2? : 


Using the 99th percentile, the maximum demand is 


X =e oa {>= 1) 
max 


ED eog04 


= bog 


TABLE If 


COMPUTATIONS IN FINDING A GOMPERTZ CURVE 


x Y Loge iteig 

5 Oa ~1.0000 

6 0.2 -0.6990 r, log ¥ = -2.0969 
jl 0.4 -0.3979 

3 0.6 =e 2ie 

9 Oey -~0.1549 penlcaai tera 37 
10 0.8 ~0.0969 

Thal 0.9 -0.0458 

2 0.9 -0.0458 r,log ¥Y = -0.0915 
13 Tae 0 


A program for the TI-59 to determine the maximum demand 
when the smallest observed X is unity, and to determine the 


parameters (k, a and b), is given in Appendix B (Program Bl). 
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Experimentation with simulated data seems to indicate that 
the Gompertz Curve procedure works well when supply is un- 
limited, but not nearly so well when supply is limited. 
EeeUoLNG THE SCHMEISER-DEUTSCH DISTRIBUTION TO ESTIMATE THE 

MAXIMUM DEMAND 

K. Pearson pointed out that a theoretical distribution must 
have at least four free parameters, in order to fit the four 
characteristics (location, dispersion, skewness, and kurtosis) 
adequately [Ref. 6]. Kottas and Lau recommends the four- 
parameter Schmeiser-Deutsch Distribution as most suitable 
(Ref. 6]. The details of the Schmeiser-Deutsch Distribution 
are given in Appendix A. 

The Schmeiser-Deutsch Distribution has a finite upper 
tail. Furthermore, if X is a random variable having the 
Schmeiser-Deutsch Distribution, there exists a closed-form 
expression for the maximum value of X. This expression is 
at+b(1-d)°, where a, b, c and d are the four free parameters. 

Using the data in the previous example, we will employ 
the Schmeiser-Deutsch Distribution to estimate the maximum 
demand. We first estimate the parameters uSing the algorithm 
in Appendix A: 

1) The mode m = (7+8)/2 = 7.5. 


Set a = 7.5. 


2) Setd (yaar y (8))7/2°= 0.5. 
3) Select the two desired quantiles which satisfy equation 


(9), namely, 
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and 


P> = 0.9, xX = ll 


Bed 4 3 


Substitute in equation (10), yielding c 
4) Substitute in equation (ll), yielding b = 52.0120. 


Then, the maximum demand is 


as 


va 
il 


a + b(l-d 
max 


14.2529 <= 14. 


This compares very well with the result obtained by using 
the Gompertz Curve which was ae = 15.3964. For these data, 
there is no substantial difference between these two alterna- 
tive methods of estimating X ax’ 

A program for the TI-59 to determine the maximum demand 
and the four parameters of the Schmeiser-Deutsch Distribution 
is given in Appendix B (Program B2). 

C. ESTIMATING MAXIMUM DEMAND WHEN THERE IS NO DATA ON 

UNSATISFIED DEMAND 

In the previous work it was implicitly assumed that the 
value of demand for each period was obtainable even if demand 
1s greater than supply. This rarely happens in the real 


world unless backorders are taken. Speaking of the "newsboy 


problem", the newsboy can collect the data in the range of 


iE, 





demand from zero to ordering quantity (Q) only. When the 
newspapers are sold out, he will not know how many more were 
in demand for that day. If he is recording the sales, his 
observed data on that day will be X = Q even if the demand 
Menactually X > Q. 

The Schmeiser-Deutsch Distribution can be easily used to 
investigate the case when data for demands greater than the 
quantity Q are not available. This is because its parameters 
can be estimated using only the mode and two selected distinct 
percentiles. 

AS an example, suppose that the newsboy orders only 10 
newspapers for each day. His sales data for a 10 day period 
Zar 6G, 9, 7, 5, 10, 10, 7, 8, 8. (Instead of demand, 
motenwewas 10, 6, 9, /, 5, 13, Il, 7, 8, 8.) The sales data 


and their cumulative relative frequency will be: 


TABLE III 


SALES DATA AND THEIR CUMULATIVE RELATIVE FREQUENCY 


x ¥(X) 
5 Okegell 
6 One 
7 0.4 
8 O26 
9 C7, 

10 1650 


Za 





Pacem aw lO) is actually Pr (X > 10), we will ignore 
this data point. 
Following the algorithm in Appendix A for the rest of the 


data points: 


II 
~ 
e 
Un 


il) Set a 


mee oet d = 0.5. 


3) Select Dia O29, = 8, oe ORE pack 9 which 


*1 2 
satisfy equation (9). 

Substitute in equation (10), yielding c = 1.5850. 

4) Substitute in equation (11), yielding b = 19.2279. 

From this, the maximum demand is X ae =— 39093 whch is 
close to ae when the newsboy has unlimited supply. 


Waen this _ one can determine Q*. Procedures to do 


this are discussed in the next chapter. 


Zee 





Tift. ESTIMATION OF OPTIMAL ORDER QUANTITY 


In the previous chapter it was shown how to estimate the 
maximum demand Cee) In this chapter, the estimate of the 
optimal order quantity (Q*) for the case of uncertainty and 


also for the case of risk will be discussed. 


A. OPTIMAL ORDER QUANTITY UNDER UNCERTAINTY 

When data in hand is not sufficient to estimate the 
possible future values of demand or when the distribution 
of demand is not known, the problem of deciding on Q falls 
into the class called decisions under uncertainty. 

It has been shown that the optimal decisions under uncer- 


maanty (O*) are {Ref. 1]: 


O 
——_ . 
Q i me ) nee , £o0r a continuous 
S Oo demand , 
and 
r C5 
as ee ar aaa * 
Q < ‘e_+e,) %max ++) I i ae Ca 


for a discrete demand, where -- 1s the unit cost of surplus 
and oe 1s the unit cost of outage. 

Since one consequence of each decision is another obser- 
vation of demand, one begins under uncertainty and uses the 
uncertainty optimal order quantity until demand data is ade- 
quate to estimate the distribution of demand and change to 


the optimum order quantity under risk. 


ae 





B. OPTIMAL ORDER QUANTITY UNDER RISK 

When one says "under risk", one means that the probability 
distribution of demand is known or, at least, some estimate 
of the probability distribution of demand is available. The 
latter case is possible using the Schmeiser-Deutsch Distri- 
bution suggested in Appendix A. 


The optimal order quantity under risk is the value of Q* 


such that 
So 
F(Q*) = —— , for a continuous demand 
Cc -- ei s s es 
Ss O GesSterivbueron, 
and 
ate) 
ko al x 
eee i) ee ane eee (ey, 
S © 


for a discrete demand distribution. 

It should be noted here that all one really needs to 
estimate the probability distribution of demand is a value 
for the [eo /(e,te.) 1 quantile. The inverse c.d.f. of the 
Schmeiser-Deutsch Distribution can easily take care of this 
problem. 

AS an example, suppose that the newsboy purchased the 
newspapers at the cost of 10¢ each and he can sell them for 
20¢ each, with leftovers having a salvage value of 5¢ each. 

Let p be the ratio [oc /(e te.) I. Therefore, p = 2/3. 

Consider the data from the previous example at the end 


of Chapter II. These data are adequate for a decision under 
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mesk. Tne ae quantile can be estimated by the inverse 
c.d.f£. of the Schmeiser-Deutsch distribution, 

ce: 
) 


a-b(d-p ne p < 4 


atb(p-d)" if pod, 


where the four parameters were found to be a = 7.5, 
Pees .22/9, c = 1.5850 and d = 0.5. 

Therefore, F7+(p) = 8.623, which is the value of Q* for 
the case of continuous demand. 

For the case of discrete demand, the optimal order quan- 
MmerEveunaer risk is the value of Q* such that F(Q*-1) < p < F(Q*), 
and then Q* = 9 in this example. 

In the next chapter we will give conclusions and applica- 


Elens from this work. 
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IV. CONCLUSIONS AND APPLICATIONS 


One of the important requirements needed to obtain the 
optimal order quantity in newsboy problems is the nature of 
the demand distribution. This thesis suggested some proce- 
dures which can be used to estimate the demand distribution 


even if data on unsatisfied demand are not available. [In 


Ss 


the case where it is hard to tell what common distribution 
the demand should belong to, the Schmeiser-Deutsch Distribu- 
tion was felt to be useful. 

It is fitting at this point to explain how the procedures 
from this and other research on the newsboy problems would 
be used. 

Suppose we are starting a sequence of single-period 
inventory ordering decisions but have no information on the 
distribution of demand. 

One way of employing these results would involve three 
phases: 

1) For the decisions at time period 1 through some period 
m, the newsboy would probably rely upon a guess at the value 
of the maximum demand (X ax)’ Considering the unit cost of 
surplus (c.) and the unit cost of outage (Co), he has two 
alternatives: 

(a) In case of C. ec Cor use Q* = X 


max 


(b) Otherwise, use the uncertainty decision rule. 


Q* ~ Ke od Wer ee "A ax 
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2) For time periods (m+l) through some later period n, 


*=s— ° , 
he would use Q [e,/(e, te.) | x with the value of xX, 


ax ax 


obtained from the Gompertz Curve or the Schmeiser-Deutsch 
Distribution, using data from the first m periods. Clearly, 
he may wish to revise this estimate as more data become 
available. 

3) After time period n, he can switch from treating the 
decision as a decision under uncertainty to a decision under 
risk, using the inverse c.d.f. of the Schmeiser-Deutsch Dis- 
tribution to estimate the value of [e /(e, +e.) 1" quantile, 
which is known to be Q* for continuous demand. 

It should be noted here that because of the data needed, 
m must provide at least three periods where he did not sell 
out, and from the work of Yong-u Sok, n should provide at 
least five periods where he did not sell out [Ref. 4]. 

From the work outlined above we have a reasonable means 
of starting a sequence of newsboy decisions when there is no 
information about the demand distribution and lost sales 
are unknown. It is hoped that the procedures suggested in 
this paper will be useful to those interested in single- 


period inventory problems and their solutions. 
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APPENDIX A 


Poole noun m EU LOCH DISTRIBUTION 


In this appendix, the Schmeiser-Deutsch Distribution 
will be briefly reviewed. For a further study [Ref. 6] and 
[Ref. 7] are suggested. 

Due to the lack of convenient methods of generating random 
values from distributions having more arbitrary shapes, the 
Schmeiser-Deutsch Distribution was developed and found suit- 
able for generating random values from many common statistical 
aroeaLoutions. 

The Schmeiser-Deutsch Distribution can take on a wide 
variety of shapes, which makes it useful for many applica- 
tions. Distributions ranging from Bernoulli trials through 
U-shaped distributions to the uniform distribution, as well 
as heavier tailed distributions, are obtainable. Also, shapes 
ranging from symmetry to maximum skewness, as measured by the 
third moment, are possible. The Schmeiser-Deutsch Distribu- 
tion can have values for the first four moments so as to 
match any of the well-known distributions, and straight- 
forward parameter determination techniques exist. 

A most favorable characteristic of the Schmeiser-Deutsch 
Distributions is the existence of simple closed-form expres- 
mmemee tor both the c.d.f. and the inverse c.d.f., being 


respectively, 


2 





d-((a-x)/b]*/° if (apace) < x <a 
a) PO D 


d+[{(x-a)/b]177° ig a < x < atb(l-a)° 


and 


a-b(d-p) © 


a+b (p-d) © if pod 


where a, b, c and d are the four parameters of the distribu- 
tion, and p is a percentile [(Ref. 7]. 

"The location and the spread of the distribution are 
determined by a and b, respectively. The shape of the dis- 
tribution, often measured in terms of skewness and tail- 


Gy. 
) in case of 


weight (or peakness), is determined by -(d-p 
Pea, Or (p-a) © Pics cmOrmeouewcd. Tnis factor 1s a trans- 
formation from the rectangular shape of the uniform distribu- 
tion to the shape of interest. The d = 0.5 yields a symmetric 
Petra oution. If c > 1, then skew is to the right for 
pmmereo and to the left for d > 0.5. (If c < l, the direction 
of skew is reversed.)" [Ref. 7] 

"The distribution possesses two characteristics which may 
adversely affect its usefulness for specific modeling pur- 


poses. The first characteristic is that the value of the 


density function evaluated at the mode, f(x = a), assumes 
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Satyeeiree Values: zero for 0 <¢c¢ <1, one force=l, 
Sreminiinity £Or c > 1. The second characteristic is the 
truncated tails of the distribution." [Ref. 7] 

"Despite the truncated tails, the distribution's versa- 
tility in assuming a wide variety of shapes makes it a 
reasonable model for a wide range of processes. This is due 
to the ability to assume the essential features of a process, 


which often makes the truncation insignificant." [Ref. 7] 


Parameter Determination [Ref. 7] 

The mode and quantiles are used in estimating parameters 
for the Schmeiser-Deutsch Distribution. The Schmeiser-Deutsch 
Distribution is useful in two common situations. The first 
is when full information 1s not available, but such quantities 
as the "most likely value" and the "minimum" and "maximum" 
values can be estimated. The second situation is when data 
is in histogram form where F(m) is the percentile of the 
mode m. The two desired quantiles, Pi = F(x,) and P> = F(x.) , 


must satisfy either of the following two expressions: 
(Ja-x,| <|[x,-al and [d-p,| <|p,-al) 


(ja-x,| <|x,-a| and |d-p, | < |p5-d]) 


1 
Although computations to determine the four parameters 

can be easily handled by computer, much attention has been 

given in this thesis to develop "manual procedures", meaning 


Manual computational capability. The TI-59 handheld calculator 
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(characterized by low price and surprising capability and 
compactness) can be considered as a "manual" tool, and a 
program for the TI-59 to determine the four parameters is 
given in Appendix B. 

An algorithm for estimating the parameters is as follows: 
(Ref. 7] 

1) Set a=m 


eee Set d = F (m) 

















By Set 
a-X d=-p 
¢ = Inj— | / In| L (10) 
3 Po 
4) Set b= |a-x,|/|d-p,|° (11) 
5) Due to the existence of the inverse c.d.f. in closed- 


form, the value of X corresponding to a percentile p of the 

distribution may be found by inserting p into equation (8). 
Furthermore, the range of X can be calculated easily 

by substituting p = 0 and p = 1 at Step 5), or obtained 

directly from equation (7) since oa = a-bd° and a = atb(1-d)°. 

At this point, one should determine whether or not the ob- 


tained distribution is an adequate model by plotting the 


fever. OL C.d.f. [Ref. 7}. 


30 





APPENDIX B 


TI-59 Programs 


This appendix presents two TI-59 programs to determine 
the maximum demand and the parameters. The first program 
is using the Gompertz Curve with the restriction that the 
smallest X is unity. The second program is using the 
Schmeiser-Deutsch Distribution with the two desired quantiles 


satisfying equation (9). 


Program BL 


i Prognam Description 


The program will estimate the maximum demand using 
the Gompertz Curve, and will also determine the parameters 


(k, a and b) in Y = ka 


Zee User Instruction 











Step Procedure Enter Press Display 
1 Enter n n [A] n 
2 Enter Yyreeer Ya, Vy Rs | vy 
ae DS Z,log 4 
Yael R/S) reed 
Yon [R75 ] Z5log Y 


ou 





Step Procedure Enter Press Display 


al a *on+l 
la (R75) Z4log 6 
3 Compute X ax (8B) Soe 
4 Review parameters 
(OpEL oma LT) 
k Kea k 
a (D] a 
b [E) b 


The maximum demand shown on the display is for the case 
when the smallest observed X is unity. Otherwise, the maxi- 
mum demand should be adjusted by adding X_. _-l, where X_. 

min min 


is the smallest observed X. 
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Program B2 
1. Program Description 


Estimate the maximum demand using the Schmeiser- 
Deutsch Distribution. Also determine the parameters (a, b, 


c and d) in 


qd - [35%) 1/s Lig a-ba® cece 4a 
F(x) = 
d + (258) 77° lp ees a+b (1-a) ° 
2. User Instruction 
step Procedure Beek Press Display 
au Enter mode m m (R7S] m 
Z Enter F(m) F (m) (R7S] F (m) 
3 Enter p, Py (R7S] P} 
+ Enter x) xy ERA xX) 
5 Enter p, P5 [R7S] P> 
6 Enter Xo Xo (R735) Xo 
7 Compute a LE] ae 
8 Review parameters 
(optional) 
a [AJ a 
b (8) b 
e [C} S 
d (DJ d 


oo 
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